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DISPLACEMENT EFFECT OF A THREE-DIMENSIONAL BOUNDARY LAYER! 


By FRANKLIN K. Moor» 


SUMMARY 


4 method is described for determining ihe “displacement 
surface” of a known three-dimensional compressible boundary- 
layer flow in terms of the mass-flow defects associated with the 
profiles of the two velocity components parallel to the surface. 
The result is a generalization of the plane flow concept of dis- 
placement thickness introduced in order to describe how a thin 
boundary layer distorts the outer nonviscous flow. 

The height of the displacement surface above the body surface 
Jor flow about a yawed infinite cylinder is shown to be equal to 
the height characterizing the mass-flow defect of the chordwise 
velocity profile. The displacement-surface height is shown to 
differ, in general, from that associated with the resultant mass- 
flow defect, even at stagnation points of the secondary flow. 
Numerical values are found for the known three-dimensional 
boundary-layer flow about a cone at a small angle of attack to 
a supersonic stream. 


INTRODUCTION 


The boundary layer established in the flow of a slightly 
viscous fluid about a body is normally considered an isolated 
region wherein the effects of viscosity predominate and 
outside of which the motion of the fluid is governed by the 
laws of nonviscous motion. For large Reynolds numbers, 
the boundary layer is assumed to be so thin that the non- 
viscous portion of the flow occurs as though there were no 
boundary layer. This assumption is strictly correct in the 
limit of infinite Reynolds number. For large but finite 
Reynolds numbers, the growth of the boundary layer causes 
the stream to be deflected away from the body surface. 

This displacement effect of the boundary layer on the 
nonviscous flow may properly be determined from the behav- 
ior of the boundary layer itself, as established either by 
experiment or by solution of the Prandtl boundary-layer 
equations for laminar flow. 

It does not follow, however, that this revised outer flow 
may properly be used in conjunction with the Prandtl 
equations to yield an improvement in the boundary-layer 
calculation. Such an improvement may be obtained only 
by use of a new set of equations that take into account the 
variation of pressure across the boundary layer. This varia- 
tion is neglected in the Prandtl equations. (See Alden’s 
iterative solution for incompressible flat-plate flow, ref. 1.) 


The customary definition of displacement thickness (ref. 2) 
is applicable to two-dimensional flow and is expressed in the 
following equation: 


[row dy=pim(h—s) (1) 


where A is some location well outside the boundary layer 
(fig. 1) at which pu= pu, and beyond which the flow may be 
considered nonviscous. (A list of symbols is provided in the 
appendix.) Under the assumption of an extremely thin 
boundary layer, h is so small that pu may be taken as the 
evaluation at the body surface (y=0) of the nonviscous flow 
obtained by neglecting the presence of the boundary layer. 
Equation (1) equates the actual mass flow near the surface 
with the mass flow which would be associated with a non- 
viscous flow that terminates at 5* rather than extending to 
the wall. Thus, the nonviscous portion of the flow behaves 
as if it occurred in the presence of a solid boundary given 
by the displacement thickness 6*(z). Equation (1) may be 


solved for 6*: 
[a 
pil 


Ordinarily, theoretical boundary-layer solutions for pu/p1%ı 
asymptotically approach 1 for large values of Reynolds 
number based on y. Therefore, displacement thickness is 
often defined as follows: 
d 
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FIGURE 1.— Plane boundary-layer flow. 


! Supersedes NACA TN 2722, “Displacement Effect of a Thres-Dimensional Boundary Layer,” by Franklin K. Moore, 1952, 
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FIGURE 2.—Three-dimensional boundary-layer flow (Cartesian 
coordinates). 


For three-dimensional boundary-layer flows, two lengths 
characterizing mass-flow defects may be defined in terms of 
the profiles of the two velocity components tangential to the 


surface (fig. 2), l 
s= fi (a e) ay | . (2) 


|, (1) e 


and it is not clear which, if either, defines a “displacement 
surface” that properly describes the extent to which the non- 
viscous flow is deflected by the boundary layer. 

Of course, it is expected that there does exist a displace- 
ment surface for such flows. The analysis that follows 
shows that such a surface may be described by a defining 
equation more fundamental than equation (1). 

The velocity and density profiles are assumed to be known 
for the three-dimensional flows under consideration. Cases 
for which this is true include: the laminar boundary layer 
on yawed infinite cylinders, treated by Prandtl (ref. 3), 
R. T. Jones (ref. 4), and Sears (ref. 5); and the laminar 
boundary layer on a cone at a small angle of attack to a 
supersonic stream (ref. 6). The displacement effect of these 
flows will be treated specifically. 

The investigation was conducted at the NACA Lewis 
laboratory in February 1952. 


THEORY 
DEFINING RELATION FOR DISPLACEMENT SURFACE 


The boundary-layer solution (assumed known) yields a 
certain distribution of velocity :(z,z) normal to the body 
surface at the outer edge h(x,2) of the boundary layer, where 
pu, and w may be taken essentially equal to p, and w 
(see fig. 2). Under the assumption that the nonviscous flow 
is altered only slightly by displacement, the most direct 
way to compute this effect would be to suppose that the 
nonviscous equations hold for y>A(x*,2) and to impose the 


following boundary condition on the normal velocity vos in 


the outer flow:, ' 
Doy=0 (x, 2) at y=h(x, 2) (3) 
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However, since the boundary conditions usually encountered 
in nonviscous flow specify am impermeable surface, it is 


convenient to recast the boundary condition given in equa- 


tion (3) in answer to the question: What impermeable 
surface y=A(z,z) would deflect a nonviscous fluid in such a 
way as to produce a normal velocity satisfying condition (3)? 
This fictitious surface may be called the displacement surface. 

Because A(z, 2) is imagined to be a stream surface, it is 


necessary to specify that at y=A the resultant velocity 


vector (u, v, w) be tangent to the surface y=A. ‘Thus, at 
y=A the ratio of » to the magnitude of q, which is defined as 


the vector (u, w), must be set equal to the slope of the surface 
y=A, measured in the direction of the vector q, or, equiva- 


lently, equal to the component in the direction of g of the 
vector grad A. In vector notation, therefore, the normal 


‘velocity » which would be produced in a nonviscous fluid at 


an impermeable surface y=A is q - grad A. The vector g 
may be obtained by evaluating the velocity vector of the 
unrevised nonviscous flow at the body surface (y=0), under 
the related assumptions that the velocity vector varies only 
slightly over distances of the order of the actual boundary- 
layer. thickness and that the revision required to take account 


of displacement is slight. Thus, at y=A (see fig. 3), 
=q . grad A 


The increment in py between A and. k is approximately 
(h-A)0v,,/Oy; again, a thin boundary layer is assumed and 
only the first term in a Taylor’s series is used. 

To the order of approximation contemplated in this 
analysis, 00,,/0y.may be obtained from the unrevised non- 
viscous flow evaluated at y=0. Thus, the fictitious im- 
permeable surface A would produce, at y=h, a normal 
velocity (see fig. 3) 

Vo= qı - grad A+(h—A) Ze) 


oy 


Introducing this result into the boundary condition given in 
equation (3) yields the defining relation for A(z, 2) 


OV, 
Codi 0—4) (FL) ge grad A (4) 
The boundary-layer solution yields vz. 


EXPRESSION OF A IN TERMS OF MASS-FLOW DEFECTS 
The displacement surface A may be related to the mass- 
flow defects (eqs. (2)), which characterize the boundary 
layer, as follows: In a Cartesian coordinate system (fig. 2), 
the equation of continuity for both the boundary-layer and 


nonviscous flow is 
Opv_ Op Opw (5) 
oy Ox 02 


Under the Prandtl boundary-layer assumptions, the Cartesian 


_ equations of motion, and hence equation (5), may be applied 


in an orthogonal curvilinear coordinate system in which the 
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suríace of the body is given by y=0, provided the de of 
curvature of the body is everywhere large as compared with 
the boundary-layer thickness. Integrating equation (5) 
across the boundary layer yields 


(P Doi)y=r = -f (+ 
= Í + E soul ni dy— 


f Opa y +2) dy 
0 or 


Set 


or, inasmuch as h is outside the boundary layer where 


pU=pU and pw= pi, 


Opt | Cert 


Pi (Voi)yar™= —h TAL 


oH) 4 (pornò) HE (ort18,) (6) 


where 6, and 6, are the mass-flow defects defined in equations 
(2). For the nonviscous flow, v,,=0 at the body surface 
(y=0), and equation (5) becomes 


Op 1 Op mt) (7) 


Mor a AS! 
P1 Oy ),_= Ox ES Oz 


Introducing equations (6) and (7) into equation (4) yields 


Š form (AH Iowa (80) 
In a Cartesian coordinate system, the displacement surface 
is related to the mass-flow defects by equation (8a). When 
cases arise for which other coordinate systems must be used, 
the following generalization in vector form of equation (8a) 
may be used: 


h 
div | pga- | (gepi [=O (8h) 


where y is the distance normal to the body surface and where 
the divergence operator involves differentiation only with 
respect to the two coordinates parallel to the body surface. 


EXAMPLES 


Plane flow.—Equation (8a) may be integrated to- yield, 
for plane flow (0/0z=0), 


A= AR 


py (9) 

where K is a constant of integration. 
The appearance of this constant means, in general, that 
the revised boundary condition on the nonviscous flow near 
the wall may be applied at any surface in the boundary-layer 
region; for example, along the wall itself. Of course, if 
there is a stagnation point on the body where u, vanishes, 


then K must be taken equal to zero and the revised boundary: 
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codos must be imposed at the location A=8,, at least 


near the stagnation point. 


dA - (52) 


Surface 
y=hlx 2) 
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` -Body surface 


FIGURE 3.—Displacement surface in three-dimensional flow (Cartesian 
coordinates). 


Stagnation point of secondary flow.—In the vertical plane 
of symmetry of the flow about a body of revolution at an 
angle of attack, the circumferential velocity component 
vanishes in the boundary layer as well as in the nonviscous 
outer region. Such stagnation of the secondary motion 
would occur in a variety of cases, in particular, wherever an 
essentially three-dimensional flow has a plane of symmetry. 
When wi is designated as the component of secondary flow 
(circumferential velocity for a body of revolution) at the 
outer edge of the boundary layer, and equation (8a) is 
written in the form 
=| pr a 5) [Hor] pra- 3) [+ 10-5) F2=0 (10) 
it is clear that A cannot equal ô+, in general, even if both w, 
and the boundary-layer ae of w vanish, unless 9w,/92 
also Vanishes. 

Flow about yawed infinite cylinder.—In the flow about a 


. yawed infinite cylinder, there is a spanwise boundary-layer 


velocity profile and an associated spanwise mass-flow defect. 
(See refs. 3 to 5.) If x is taken in the chordwise direction 
(fig. 4), the entire flow depends only on +. ‘Thus, derivatives 
with respect to the spanwise coordinate z vanish; and, hence, 
from equation (8a), the plane-flow result (eq. (9)) applies. 
Accordingly, the spanwise mass-flow defect represented by 
ö, does not enter into the determination of the displacement 
surface. 


FIGURE 4.—Coordinate system for yawed infinite cylinder. 
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FiqGurE 5.—Coordinate system for cone at angle of attack. 


Supersonic flow about cone at small angle of attack.— 
When the coordinate system shown in figure 5 is used, 
equation (Sb) becomes 


sin 9 | orna ae) es ES (A—8,) |=0 (11) 


where i 
5,= 1-27) d 
f f, ( piui y "a 
A $ aD 
j= i ) a 
mE Í, ( piwi y 


Because the outer flow is conical, pı, Wt, and w, are functions 
only of »; and equation (11) may be written 


(12) 


pr sin 03 26-52 Hi| pman |=0 (18) 


In the case of conical outer flow, the associated boundary- 
layer profiles show similarity of the Blasius type in meridional 
planes (see ref. 7 or 8). Thus, in a meridional plane, A, ôz, 
and ô are proportional to yz. Incorporation of this infor- 
mation into equation (13) gives 


; Ò 
5 pit; sin AI E (459) |=0 


For a cone at small angle of attack a, u, is nearly equal to , 


u, the velocity on the cone surface at zero angle of attack. 
The quantities pı, A, ö,, and 3, vary only slightly with angle 
of attack, whereas 

Wy = aA sin. p 


where A, depends only on the cone vertex angle and the flow 
Mach number and is defined in reference 6 as follows: 


(15) 


The quantities z, x, and Y are in the notation of reference 9, 
wherein they are tabulated. E 


(14) 
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To the first order in a, the substitution of equation (15) 
into equation (14) yields the result 


Lin As (ö.—5.) COS Y (16; 


Amita sin 6 


The analysis of reference 6 yields the values of ö, and 6). 
Clearly; A differs from 6, in the plane of symmetry p=0, r, 


where the circumferential velocity w vanishes. It might, 
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Fısur& 6.—Proportional rate of increase of displacement thickness 
with angle of attack on cone. 
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2-2 
6 shows the proportional rate of increase of displacement 
thickness A with angle of attack in the plane of symmetry 
y=r at zero angle of attack. The corresponding rate of 
increase of the mass-flow defect 5, is shown for comparison. 
These curves are obtained from equation (16) and the 
results of reference 6. The sine of the semivertex angle is 
introduced as a factor primarily to permit presentation of 
the case 6=0 as the limit of an indeterminate form. For a 
stream Mach number of 2, the change in displacement thick- 
ness appears to be of the order of 50 percent larger than the 
change in the meridional mass-flow defect. 


however, be noted that in the plane p= , A=6,. Figure 
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CONCLUDING REMARKS 


The foregoing analysis deals only with the displacement 
effect of a known boundary layer on the nonviscous outer 
flow, and hence can be applied only if the boundary-layer 
behavior has been determined either theoretically or experi- 
mentally, The latter approach might possibly find appli- 
cation in the correction of nozzle contours for boundary-layer 
development. In this connection, perhaps, it should be 
noted that the analysis is not restricted to laminar flows. 


Lewis FLIGET PROPULSION LABORATORY 
NATIONAL ADYISORY COMMITTEB FOR ÁERONAUTICS 
CLEVELAND, Onto, March 6, 1952 


APPENDIX—SYMBOLS 
The following symbols are used in this report: l : 
Ay function of cone angle and Mach number | z coordinates in body surface - 
(eq. (15)) Y coordinate normal to surface 
h height above body surface at which p, u, | « angle of attack 
W= pi, UA, Wi, and beyond which nonviscous | A ‚height above body of displacement surface - 
equations apply (eq. (3)) 5* displacement thickness in plane flow 
M, stream Mach number bs length characterizing mass-flow defect of 
q velocity vector composed of components paral- u-profile (eqs. (2) and (12)) 


lel to body surface u, w 


bs length characterizing mass-flow defect of 


U velocity component in x-direction w-profile (eq. (12)) 
U meridional velocity component at surface of | 3, length characterizing mass-flow defect of 
cone at zero angle of attack w-profile (eq. (12)) Ä 
v velocity component in y-direction 8 semivertex angle of cone 
Vor boundary-layer solution for velocity normal to | p density 
surface, evaluated at outer edge of boundary | » angular coordinate (fig. 5) 
layer . 
Vos nonviscous solution for velocity normal to | Subscript 1 denotes evaluation of nonviscous flow at body 
surface, evaluated near the surface surface, taken equivalent to conditions at outer edge of 
w velocity component in z-direction boundary layer of infinitesimal thickness. 
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